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ABSTRACT 

Observational studies have revealed a "downsizing" trend in black hole (BH) 
growth: the number densities of luminous AGN peak at higher redshifts than those 
of faint AGN. This would seem to imply that massive black holes formed before low 
mass black holes, in apparent contradiction to hierarchical clustering scenarios. We 
investigate whether this observed "downsizing" in BH growth is reproduced in a semi- 
analytic model for the formation and evolution of galaxies and black holes, set within 
the hierarchical paradigm for structure formation (Somerville et al. 2008; S08). In this 
model, black holes evolve from light seeds (~ lOOAf ) and their growth is merger- 
driven. The original S08 model (baseline model) reproduces the number density of 
AGN at intermediate redshifts and luminosities, but underproduces luminous AGN at 
very high redshift (z > 3) and overproduces them at low redshift (z < 1). In addition, 
the baseline model underproduces low- luminosity AGN at low redshift (z < 1). In or- 
der to solve these problems, we consider several modifications to the physical processes 
in the model: (1) a 'heavy' black hole seeding scenario (2) a sub-Eddington accretion 
rate ceiling that depends on the cold gas fraction, and (3) an additional black hole 
accretion mode due to disk instabilities. With these three modifications, the models 
can explain the observed downsizing, successfully reproduce the bolometric AGN lu- 
minosity function and simultaneously reproduce galaxy and black hole properties in 
the local Universe. We also perform a comparison with the observed soft and hard 
X-ray luminosity functions of AGN, including an empirical correction for torus-level 
obscuration, and reach similar conclusions. Our best-fit model suggests a scenario in 
which disk instabilities are the main driver for moderately luminous Seyfert galaxies 
at low redshift, while major mergers are the main trigger for luminous AGN. 
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1 INTRODUCTION 

It has been known since early optical quasar surveys 
that the co-moving number density of luminous quasars 
has a pronounced peak at a redshift around z = 2- 
2.5 dSchmidt fc Greenlll983l;lBovle et al.lll98Sl;lHewett et al 



1994 : iBovle et al.ll200d : I Warren et all 1 19941 : Ischmidt et al 
1995), with a fairly steep decline at higher and lower 



redshift. More recently, it was discovered that very mas- 
sive B H (10 9 M(T)) appe a r to e xist already at z ~ 6 and 
higher jFan et all l200d . l200ll : iMortlock et~afl 120111 ). but 
they are extremely rare, in accordance with this trend. 
Recent progress in detecting faint and obscured AGN 
has been achieved by analysing data from X-ray surveys 
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(XMM-Newton, Chandra, ROSAT, ASCA, e.g. 


Mivaii ct al. 


200d: lLa Franca et al.l 20021: Cowie et al. 2002 


; Fiore et al. 


2003: Barger et al. 2003: lUeda et al.l 2003; Hasineer et al. 


20051: iBareer & Cowid 2005; ISazonov & Revnivtsev 


2004; 


Nandra et all 120051; Ebrero et al. 20091; Aird et al. 


2010; 


Fiore et al.ll2012f). All of these studies in the hard and soft X- 



ray range find that the cosmic evolution of AGN is strongly 
dependent on the AGN luminosity: the number density of 
successively less luminous AGN peaks at lower redshifts, 
with the lowest luminosity AGN showing a basically con- 
stant number density. Making the simplified assumption 
that AGN luminosity is proportional to BH mass (as we 
would expect if black holes are accreting at the Eddington 
rate, L oc M.) would imply that very massive black holes 
seem to be already in place at very early times, whereas 
less massive black holes seem to evolve predominantly at 
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lower redshifts. This behavior is called 'downsizing' or 'anti- 
hierarchical' growth of b lack holes. The down s izing trend is 
also seen in the optical llCristiani et al]|2004 ICroom et "aD 
2004lFan et al.ll2004lHunt et al]|2004lRichards et alJIgOOfj ; 
Wolf et al.ll200a) and the NIR (e.g. iMatute et al.ll2006ft .~ 



On the face of it, this observational result seems to be 
in conflict with the expectations within the currently fa- 
vored hierarchical structure formation paradigm, such as 
those based on the Cold Dark Matter (CDM) model l|Peeblesl 
1 19651 ; IWhite fe Reesl Il978l ; iBlumenthal et alj 1 1985ft . In this 
framework, low mass halos form first and more massive halos 
grow over time via subsequent merging and smooth accre- 
tion. However it is now well known that the evolutionary his- 
tory of observable galaxies also follows an "anti-hierarchical" 
or downsizing behavior, with several independent obser- 
vational indicators suggesting that more massive galaxies 
formed their stars and had their star formation "quenched" 
earlier than low-mass galaxies, which continue forming stars 
to the present day (an o verview of these observations is given 
in lFontanot et aij|200g| ). 

Present-day spheroidal gala xies host supermassive 
black holes at thei r centers (|Magorrian et aL 1998; 
iGenzel fe Eckartl 1 1999ft and strong correlations have been 
found bet ween black hole masses and properties of the i r host 
galaxies ilFerrarese fe Merritd |2000|; [Gebhardt et alJ l2000l ; 
iTremaine et all |2002| ; lHaring fe Rixl hoo4 iGrahaml \20m . 
This can be interpreted as evidence for co-evolution between 
the host galaxies and their black holes, but some observa- 
tions of BH growth and SF in individual objects appear 
to contra dict the picture of sim ple one-to-one co-evolution 
over time l|Mullanev et al.ll2012h . Thus the details of BH and 
galaxy co-evolution remain unclear. 

During their lifetime, black holes are assumed to un- 
dergo several episodes of significant gas accretion, dur- 
ing which this accretion p owers luminous quasars or ac- 
tive galactic nucle i (AGN) (|Salpeterlll96l ; IZel'Dovich|[l964 
iLvnden-Bellll 1969ft . By estimating the total energy radiated 
by AGN over their whole lifetime, it can be shown that 
nearly all the mass seen in dormant black holes today can 
be accu mulated durin g the periods of observed bright AGN 
activity |Soltanlll982T ). This implies that there is not a great 
deal of room for "dark" or obscured accretion. 

A large number of works have explored the predic- 
tions of the ACDM model for the formation and evolu- 
tion of supermassive black holes and AGN, with vary- 
ing levels of complexity. Several works made predictions 
based on nearly purely analyt i c models dEfstathiou fe Reesl 
1 19881 ; lHaehnelt fe Reesl Il993l; lHaiman fe Loebl Il998ft and 
on semi-empirical models (jShankar et alj|2009t I2010L 120121 ') 
and a large number of studies have been published 
based on semi-analytic models of galaxy formation within 
which mechanisms for black hol e formation and evo- 

lution have been included (e.g. Kauffmann fe Haehnelt 
-——1 r-r- : n U^ol — — 1 J L T 



2000; Volontcri ct al. 2003; Granato et al. 2004; Mcnci ct al 
2004 iBromlev et alj|2004 ICrotoiJl2006l: iBower et al.ll2006l : 



Marulli et alj |2008l ; ISomerville et all 120081 ). Recently, nu 
merical hydrodynamic simulations have also included black 
hole growth and AGN feedback using "sub-grid" recipes 
ilSpringel et al.ll2005l; iHopkins et~aH 120061; Ibi Matteo et al ' 



Hopk ins et al] (|2008l ) presented predictions based on "semi- 
empirical" models in which galaxy properties were taken 
from observations and the relationship between galaxy prop- 
erties and AGN luminosity was based on the results of a 
large suite of hydrodynamic merger simulations. 

These models differ in many of the details of how BH 
formation and growth are implemented, but there seems to 
be a broad consensus on several points. First, self-regulated 
black h ole growth, perhaps v ia radiation pressure driven 
wind s (iDi Matteo et al] 120051 ; iKind 120051 ; iMurrav et ail 
120051; iRobertson et al.ll2006ft is a widely adopted means of 
obtaining the observed tight relationship between BH mass 
and spheroid mass (although some models simply assume 
this relationship without invoking a physical mechanism). 
Second, one must invoke a physical mechanism that can 
feed large amounts of gas onto the central BH within a 
short time. Galaxy-galaxy mergers are a popular (though 
not universally adopted) way to remove angular momen- 
tum from the gas and efficiently drive it to the center 
of the galaxy. Third, some mechanism must be adopted 
that reduces or stops gas cooling and accretion in mas- 
sive dark matter halos. Many models assume that low levels 
of accretion onto SMBH can produce radio jets that heat 
the surrounding hot h alo i|Crotonll20o"6l ; IBower et al.ll2006l ; 
ISomerville et al]|2008ft thereby inhibiting cooling flows, pre- 
venting over-massive galaxies from forming and quenching 
star formation in massive galaxies, leading to more obser- 
vationally co nsistent color-magn itude or SFR-stellar mass 
relationships (|Kimm et al.ll2009ft . 

However, several aspects of the physics of black 
hole growth and formation remain poorly understood. 
First, there are active ongoing debates about when and 
how seed black holes form: either via a direct col- 
lapse of cold g as clouds leading t o massive seeds of 
10 4 -10 5 M W dLoeb fe Rasid 1 1994 iKoushiappas et~aH 



20051 ; Robertson et all 120061 ; iLi et alj 120071; ISiiacki et al. 
20071; iJohansson et alj 120091; iMcCarthv et alj |2010| l201ll : 



Degraf et alj|201ll ; lDi Matteo et alj|2012l ; IChoi et al] |2012) 



|2004 IVolonteri fe Starkl l20liF iBellovarv et alJ l201lft or 

via stellar remnants from Pop III s tars, resulting in low 
mass (~ lOOMp)) blac k hole seeds (|Madau fe Reesl l200ll ; 
iHeger fe Wooslevll2002ft . An alternative possibility is direct 
seed formatio n in a merger even t as seen in the numerical 
simulations of iMaver et alj (|2010ft . While the seeding mech- 
anisms do not alter the AGN and black hole population at 
low redshift (as gas accretion during the evolution exceeds 
the seed black hole masses by many orders of magnitude), 
at high redshift the choice of the seeding model strongly 
influences the black hole formation and is highly relevant 
for understanding the observed population of luminous high 
redshift quasars. 

Another matter of vigorous debate is the process or 
processes that trigger and regulate accretion onto the cen- 
tral SMBH. As noted above, the most luminous observed 
quasars require accretion rates such that ~ 10 8 -10 9 M© of 
gas must be funneled onto the black hole over a timescale 
of < 10 8 ' 5 yr, i.e. nearly the whole gas content of a good- 
sized galaxy must be fed onto the b lack hole in roughly 
a dynamical time (see discussion in IHopkins" et all 120081 . 
and references therein). Mergers appear to be a physically 
well-motivated candidate for producing this dramatic effect, 
and semi-empirical calculations have shown that there is 
a statistical consistency between the observed merger rate 
and the observed AGN dut y cycle and luminosity function 
l|Hopkins et alj 120061 . | 20081 ) . suggesting that it is at least 
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possible that these processes are causally linked. However, 
the observational situation remains murky. Observational 
studies have repeatedly failed to find evidence for a statisti- 
cally significant enhancement of merger-related signatures, 
such as close pairs or morphological di s turbance, in AGN 
hosts up to z ~ 1 (ICisternas et alJl20ld: iGeorgakakis et all 



20091; IPierce et al.ll2007l ; iGrogin et alj|2005l ; iLi et al.l 120081: 



Ellison et al.l 120081 ). But a recent study by lEllison et al.l 
1 20111 ) does find a significant enhancement of AGN activity 
in close pairs, and discusses reasons that previous studies 
based on similar data have not found a signal. The absence 
of enhanced merger signatures in morphological studies has 
recently been shown to persist in X-ray selected AGN up to 
z ~ 2.5 l|Schawinski et alj|201lt iKocevski et al.ll2012f ). Re- 
cent work suggests, though, that the fraction of hosts with 
morphological disturbances may be higher in obscured AGN, 
many of which are missing from X-ray selected surveys (S. 
Juneau, private communication). 

The goal of this work is to explore the the interplay 
of different physical processes that determine the masses of 
seed black holes, the triggering of AGN activity and the ef- 
ficiency of gas accretion during the active phases of black 
holes, with the aim of understanding the physical origin of 
the observed downsizing trend in black hole growth. We fol- 
low a semi-analytic approach, which has been shown to suc- 
cessful ly reproduce many ob s erved galaxy populatio n prop- 
erties (|Somerville et aUl2008l . lSomerville et aLlfeOHl ) and in- 
cludes a model for the merger triggered formation and evo- 
lution of black holes (see Section [3] for details). In this paper 
we present three major modifications to the baseline model 
published in S08: 

• gas fraction dependent Eddington ratios (accretion ef- 
ficiency) 

• AGN activity triggered by disk instabilities 

• 'heavy' black hole seeds 

to the existing model. We find that with these three modifi- 
cations, our model does reproduce the observed downsizing 
behavior. 

In Section [2] we briefly summarize the results from pre- 
vious studies. Section [3] gives a brief overview of the semi- 
analytic model used in our study and describes the different 
modifications for black hole growth we are considering. In 
Sections [4] and [5] we present some properties of present-day 
and high-redshift galaxies and their black holes, which are 
compared to observations. The AGN number density evolu- 
tion is studied in Section [6] considering the influence of the 
different modifications concerning black hole growth. Sec- 
tion[7]presents a comparison of the evolution of the observed 
bolometric and hard and soft X-ray luminosity function with 
the model output. In Sections [5] and [5] we discuss the evolu- 
tion of the Eddington ratio distributions and the evolution 
of the black hole- AGN luminosity plane. Finally, in Section 
1101 we summarize and discuss our main results. 



2 PREVIOUS STUDIES 

A number of studies have investigated the observed 
anti- hierarchical tre nd of BH a c tivity using the 'Gal- 
form' (Durham) llBower et al.l 120061 ). the 'Munich' 



IIDe Lucia fc Blaizodl2007T l or the 'MORGANA' semi-analytic 



models (|Monaco fc FontanotJ 120051 ; iFontanot et~akl l2006h . 
The first two models are applied to the dark matter merger 
trees of the Millennium simulation, while the latter uses 
merger trees from the PlNOCCHlO method. All models dis- 
tinguish between black hole accretion in the bright quasar 
mode and the low- Eddington ratio radio-mode. In the 'Gal- 
form' model, the quasar mode is triggered by merger events 
and disk instabilities, but the mass growth of black holes is 
dominated by accretion due to disk instabilities. In the 'Mu- 
nich' model the quasar mode is assumed to be triggered only 
by merger events. The 'Morgana' model assumes that any 
low angular momentum gas within the "bulge" gas reser- 
voir is available to accrete onto the central BH. Gas may 
be transferred to the "low- J" reservoir either by mergers 
or disk instabilities. In this model, black holes grow at a 
rate determined by the viscosity of the accretion disk and 
the mass of gas in the low-J reservoir. In the 'Galform'- 
model the black hole accretion rate is proportional to the 
star formation rate during a starburst (triggered either by a 
merger or disk instability), whereas in the 'MuNlCH'-model 
the accretion rate is dependent on the cold gas content in 
the galaxy, the galaxy circular velocity, and the mass ratio 
of the triggering merger. While these recipes differ in detail, 
they are both proportional to the gas content and rougly 
inversely proportional to (1 + Vc -2 ), where V c is the circular 
velocity of the bulge component in the case of 'Galform' 
and the halo virial velocity in the 'MuNlCH'-model, and both 
lead to a black hole-bulge mass relation at z = that is con- 
sistent with the observed one. However, it appears that the 
history of black hole accretion is significantly different in the 
two models. In the ' MORGANA'-model, black holes grow at 
a rate determined by the viscosity of the accretion disk and 
the mass of gas in the low-J reservoir. 

IFontanot et al.l l|2006l ) claim that their 'Morgana'- 
model can reproduce downsizing, which in their model is 
caused by stellar kinetic feedback that arises in star-forming 
bulges leading to a removal of cold gas in small elliptical 
galaxies (reduction of the number of faint AGN at high red- 
shift). To obtain a good match to the number density of 
bright quasars they require quasar-triggered galactic winds, 
which self- limit the accretion onto black holes. 

iMalbon et all (|2007l ) found, using the 'Galform' semi- 
analytic code, that the direct accretion of cold gas during 
starbursts is an important growth mechanism for lower mass 
black holes and and for all black holes at high redshift. The 
assembly of pre-existing black hole mass into larger units 
via merging dominates the growth of more massive black 
holes at low redshift. Therefore, they claim that as redshift 
decreases, progressively less massive black holes have the 
highest growth rates, in agreement with downsizing. Their 
model output reproduces the evolution of the optical lu- 
minosity function of quasars, however, they do not show 
a quantitative comparison for the X-ray and/or bolometric 
AGN luminosity. 

iFanidakis et alj (|2010l ) have constructed a model based 
on the 'Galform '-model framework, but with different 
recipes for black hole growth and AGN feedback. They 
present a quantitative comparison of their model output 
to the observed quasar luminosity function at different red- 
shifts. The previous 'GALFORM'-based SAMs associated ra- 
diatively efficient BH accretion with "cold mode" accretion 
(merger or disc-instability triggered) and radiatively ineffi- 
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cient accretion with "hot mode" accretion ( accretion from a 
quasi- hydrostatic hot gas halo). In contrast. [Fanidaki s et al.l 
(2010) associate accretion at low accretion rates (less than 
about one percent of the Eddington rate), regardless of 
its origin, with radiatively inefficient advection dominated 
accretion flows (ADAF), and assume that accretion at 
higher rates wi ll be r adiatively efficient. At high redshift, 
iFanidakis et al.l (|2010T ) do not Eddington-limit the accre- 
tion rates, but allow super-Eddington accretion, which is 
responsible for very luminous AGN at high redshifts. The 
high number density of low luminosity AGN at low redshift 
can be explained by the luminosity produced via the ADAF 
mode. They attribute the observed downsizing trend mainly 
to d ust obscuration of lo w luminosity AGN at high redshift. 

iMarulli et al.l (2008) investigated different parameteri- 
zations of the quasar light curves in the 'MuNlCH'-model, 
and found as expected, that the lightcurve parameteri- 
zation has a large effect on the number density of faint 
AGN as a function of redshift. They found the best re- 
sults with a lightcurve model that includes an Eddington 
growth ph ase followed by a pow er-law decline phase as sug- 
gested by iHopkins et all l|2006h . However, for all adopted 
lightcurve models they found that the previously published 
Munich model underpredicts the number density of lumi- 
nous quas ars at high redshift (z > 1). Finally, in a follow-up 
study by iBonoli et all (|2009l ). the BH accretion efficiency 
was assumed to be a function of redshift as well as gas con- 
tent and merger mass ratio. They obtained improved results, 
but still were not able to reproduce the luminosity function 
of observed AGN over the full range in redshift and lumi- 
nosity. 



3 THE SEMI-ANALYTIC MODEL 

The semi-analytic model used in this study is presented in 
S08 and we refer the reader to this paper for details. The 
galaxy formation model is based on dark matter merger 
trees generated by the extended Press-Schechter formal- 
ism. The evolution of baryons within these dark matter 
halos is modeled using prescriptions for gas cooling, re- 
ionization, star formation, supernova feedback, metal evo- 
lution, black hole growth and AGN feedback. Here we focus 
on the mechanism describing the formation and evolution of 
black holes. Each top-level dark matter halo is seeded with 
a 100 Mq black hole in its center, which can grow by two 
mechanisms: through cold gas accretion during the 'bright' 
quasar mode and through accretion of gas from the hot 
halo via a cooling flow during the low-Eddington ratio and 
radiatively inefficient radio mode. The quasar mode is as- 
sumed to produce momentum-driven winds, which are mod- 
eled using the analytic scaling derived an d calibrated from 



binary hydrodynamic merger simulatio ns (Robertson ct al 



20061; ICox et al.ll2006al;lRobertson et al.ll2006l ; lHopkins et al 
2007bl ; iRobertson et al.ll2006l ) 



The radio mode can only occur when a hot quasi- 
hydrostatic halo is present, which is assumed to be the case 
when the cooling shock is predicted to be within the virial ra- 
dius (r coo i < r v i r ). During t his phase, bl ack holes are fuelled 
by Bondi-Hoyle-accretion {Bondi 195 2]), with t he is other- 
mal cooling-flow solution from Nulsen fc Fabian] (|2000h . The 
growth in the radio mode is also associated with an efficient 



production of radio jets that results in an energy injection 
into the intracluster medium (ICM). Therefore, we assume 
that the energy arising from the accretion onto the black 
hole couples to and heats the gas in the surrounding hot 
halo (Radio mode feedback). 

3.1 Standard accretion model 

The quasar phase is triggered by galaxy merger events with 
a mass ratio of [i > 0.1, where [i is the mass ratio of 
the baryonic components and the dark matter within the 
central part of the galaxy (see S08 for the precise defini- 
tion) . The lower limit is motivated by binary hydrodynamic 
merger simulations. Whenever the two progenitor galaxies 
merge, their black holes are assumed to also merge and 
form a single black hole whose mass is the sum of the pro- 
genitor BH's masses. The model for gas accretion onto the 
black hole is motivated by the analysis of gas inflow rates 
onto t he nuclear regions from idealized disk merger simula- 
tions llSpringel et al.ll2005l; 



l2006bl ; IHopkins et all 12*006 



Robertson et al.ll2006l ; TCox et all 



2007aj). During the merger, the 



BH is assumed to grow rapidly with accretion rates near the 
Eddington limit. This rapid accretion continues until the 
energy being deposited into the ISM in the central region 
of the galaxy is sufficient to significantly offset and eventu- 
ally halt accretion via a pressure-driven outflow. During this 
"blow-out" phase, the accretion rate declines gradually until 
the nuclear fuel is exhausted. 

Based on the merger simulations, the final black hole 
mass Af.^flnai at the end of the blow-out phase is assumed 
to be related to the mass of the spheroidal component after 
the merger: 



Here, M spn is the final spheroid mass after the merger, 
/bh, final is an adjustable parameter and T(z) describes the 
evolution of the blac k hole-bulge mass relation with time 
(Hopkins et al.i r2006X Following the merger simulations, a 
Gaussian distributed scatter with a value of a. — 0.3 dex is 
additionally applied to the accreted gas mass, representing 
stochasticity due to e.g. the properties of the orbit. When 
the black hole mass has reached its final mass value, the 
quasar mode is switched off. During the quasar phase, the 
light curve models describe two different growth regimes: 
an Eddington-limited and a power-law decline phase of ac- 
cretion. In the first regime, the black hole accretes at the 
Eddington limit until it reaches a critical black hole mass 

Af.crit: 



M., oli t = /BH.crit X 1.07 (Mm , final J 



(2) 



Here, the parameter /BH,crit is set according to the merger 
simulations, and determines how much of the black hole 
growth occurs in the Eddington-limited versus power-law 
decline phase. The growth of the black hole during the first 
regime can be modeled by an exponential increase of mass: 



M., new (t) = M. exp /edd" 



1 - e 



t 



t. 



alp 



(3) 



where e = 0.1 is the efficiency of the conversion of rest mass 
to energy, and t sa i p ~ 0.45 Gyr is the Salpeter timescale. No 
strong observational constraints are available for e and if or 
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how it evolves with redshift. However, some observations a t 
z = indicate that 0.04 < e < 0.16 |Marconi et al.ll2004h . 
For simplicity we take a constant mean value of e = 0.1 
at all redshifts, which is a very standard assumption. The 
parameter /Edd in Eq. [3] is the maximum accretion rate, 
which is defined by the ratio of bolometric luminosity to the 
Eddington luminosity: 



/Edd — ^bol/^Edd- 



(4) 



The Eddington luminosity I/Edd (assuming a hydrostatic 
equilibrium between the inward gravitational force and the 
outward radiation pressure) is given by: 

L E dd = — = 1.4 x 10 ( ] erg/s, (5) 



0°T 



Vio 8 M Q 



where <tt is the Thomson cross section for an electron and 
m v the mass of a proton. Combining eq. [3] [4] and [5] and the 
relation L\, \ = e/(l — e)Mc 2 , the corresponding accretion 
rate in the first regime can be calculated by: 

M.,i(t) = 1.26 x 10 38 erg/s — %M. lMW («). (6) 

Note that in the baseline model, the maximum accretion 
rate is assumed to equal the Eddington rate, i.e. /Edd = 1. 

Once it exceeds the critical mass M, lCr it in Eq. [2] the 
black hole enters the second regime, the 'blow-out' phase, 
which is described by a power-law decline in the accretion 
rate. Fitting the lig ht curves in merger simulations from 
iHopkins et "all (|2006T l gives the following parametrization for 
M.,u- 



M.,j/(t) 



M., pe 



i + it/t Q y 



(7) 



where £q oc i sa i p is the e-folding time, M, tPca ^ is the peak 
accretion rate (given by /Edd times the Eddington accretion 
rate) and /3 is a parameterized function of the peak accretion 
rate. In the case that the initial black hole mass is already 
larger than the calculated critical mass, the black hole is not 
allowed to accrete at the Eddington rate at all and goes im- 
mediately into the blow-out phase. If the initial black hole is 
even larger than the calculated final mass, no quasar phase 
occurs at all. Note that for calculating the bolometric lumi- 
nosity only the accretion rates during the quasar phases are 
taken into account (thus ignoring the contribution from the 
radio mode accretion): 



M.,qso c 



(8) 



where M. : qso = M.j in regime I and M.,qso = M.j in 
regime II. 



3.2 Sub-Eddington limit for the maximum 
accretion rate 

Observational studies show that the peak in the Ed- 
dington ratio distributions of QSOs is not constant with 
time; instead, it is found to be dep endent on red- 
shift as well as on black hole mass dPadovanil 1 19891; 
Vestergaard 20031: Shankar et al. [2004 ; Kollmeier et al.l 



20061: INetzer fe Trakhtenbrod 120071 : Iffickox et all 1200a 
Schulze fc W isotzki 2010). In particular at low redshifts 
Z < 1, it has been claimed that there is a sub-Eddington 



limit for black hole accretion, which is dependent on black 
hole mass and redshift (e.g. INetzer fc Trakhtenbrod 120071 ; 
ISteinhardt fc Elvidlioioh . The underlying reason for such 
a sub-Eddington limit is not obvious but might be re- 
lated to the cold ga s content of the galaxy. For example, 
IHopkins et all (2008) found, in their semi-empirical models, 
that allowing 'dry' (gas-poor) mergers to trigger quasar ac- 
tivity would overproduce luminous quasars at low redshift 
(z < 1). To explore this effect we introduce a limit for the 
Eddington ratio at z ^ 1, which is dependent on the cold 
gas fraction / co id = M co id/(M co id + M s teliar) of the merged 
galaxy after the merger. For / co id > 0.3, we still allow the 
black hole to accrete up to the Eddington-rate, while for 
/cold < 0.3 we assume a simple, linearly decreasing function 
for the maximum accretion rate: 



/Edd(/cold) = 3.3 x / co ld + 0.001 



(9) 



This lowers the peak accretion rate, and therefore also cor- 
respondingly decreases the accretion rate in the power-law 
decline part of the light curve. Note that the assumption of 
a limited accretion rate is also supported by t he results of 
semi-empirical models bv lShankar et afl (|20 1 lT ) . which also 
favor a decreasing Eddington-ratio with time and a radiative 
efficiency increasing with black hole mass. 



3.3 Disk instabilities 

Various observational studies suggest that moderately lumi- 
nous AGN are typica l ly not major-merger driven, at z < 1 



JCisternas et al. 2010l; iGeorgakakis et al.ll2009l ; [Pierce et all 
120071 ; iGrogin et al.l l2005l; ISalucci et aLl ll999T)7and interest- 
ingly also at z « 2 dKocevski et al.l l2012; Ro sario et al.ll201ll ; 
ISilverman et alj l201ll )~ as they do not find more morpho- 
logical distortions for AGN host galaxies than for inactive 
galaxies. This suggests that moderately luminous AGN may 
undergo a 'main sequence' secular growth, e.g. their nuclear 
activity might be additionally driv en by disk instabilities . 
Here we use a statistic proposed bv lEfstathiou et all (Il982f ) 
to quantify disk stability, based on numerical N-body simu- 
lations. They find that the disk becomes unstable if the ratio 
of dark matter mass to disk mass becomes smaller than a 
critical value, and give the following parameterization for 
the onset of disk instabilities: 



Mdi 



disk.crit 



G e 



(10) 



where Mdisk.crit is the critical disk mass, above which the 
disk is assumed to become unstable, w max is the maximum 
circular velocity, -Rdisk the exponential disk length and e 
the stability parameter. We use a slightly smal ler value 
(e = 0.75) than proposed in lEfstathiou et all l|l982h . as it re- 
sults in a better quantitative match with the observed AGN 
luminosity function (see section [6]). However, we find that 
the number density of AGN is not very sensitive to the pre- 
cise value of the stability parameter. A smaller stability pa- 
rameter means that the critical disk mass becomes slightly 
larger and thus, on average disks become unstable at a later 
time, leading to a minor decrease in the number density of 
instability-driven AGN. Moreover, we have seen that recent 
simulations of isolated disk galaxies tend to indicate an even 
smaller stability parameter of e ~ 0.6. Furthermore, in our 
model it is assumed that whenever the current disk mass 
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(consisting of both stars and cold gas) exceeds the critical 
disk mass, the bulge component is enlarged by the differ- 
ence of the 'excess' stellar mass AAf d i sk = Afdisk — Mdisk.crit 
so that the disk becomes stable again. We assume that a 
certain amount of cold gas (proportional to the excess disk 
mass) is additionally accreted onto the black hole triggering 
an active phase: 



M dis 



x AM dis 



(11) 



Here, we adopt /BH,disk = 10~ 3 , motivated by the local 
black hole-bulge mass relation. For the accretion process we 
assume a constant Eddington ratio of / cd d = 0.01 with a 
Gaussian distributed scatter of 0.02 dex (motivated by ob- 
servations, D. Alexander private communication). The black 
hole accretes as long as there is gas fuel from disk instabili- 
ties left (Mdiskfuci > 0) and the accretion rate is calculated 
by: 



M., disk = 1.26 x 10 38 erg/s 



e /e. 



-M. 



(12) 



For the total bolometric luminosity, the bolometric lumi- 
nosities from the merger driven quasar phase (equation [5J 
and from any disk instablities are summed up: 



+ M.,disk I c 



(13) 



In contrast, the studies using t he Munich model 
ijMarulli et al.ll2008l : iBonoli et alJl2009l ) do not consider disk 
instabilities for calculating the AGN bolometric luminosity. 
In the Galform model, a similar approach is used to es- 
timate when a disk should become unstable, but it is then 
assumed that all gas and stars in the disk are transferred 
to the bulge component, leading to a much more dramatic 
effect. As a result, in their model, disk instabilities are found 
to be the major physica l process responsible for black hole 
grow th at all redshifts (|Bower et al.l l200rj : iFanidakis et alj 
l2O10h . 

Our simple model is based on simulations of isolated 
disk galaxies, which develop secular internal instabilities. 
These secular instabilities are not expected to be associ- 
ated with large nuclear inflows nor with dramatic mor- 
phological or dynamical transformation. However, in a cos- 
mological context, the expected rapid inflows can lead to 
more violent "stream fed" disk instabilities, particularly at 
high redshift, which may drive BH feeding at high rates, 
and more dramatic morphological /dynamical t ransforma- 
tion (|Bournaud et al.ll201ll : [Ceverino et al.ll20ld ). However, 
it is not known how common such violent instabilities might 
be in a cosmological context, nor how to model their effects 
within a semi-analytic model. This is an important topic for 
future work, but for the moment we restrict ourselves to the 
more mild secular instabilities. The reader should keep in 
mind, however, that this may represent a minimal predic- 
tion for the impact of internal instabilities on AGN feeding 
and spheroid formation. 

3.4 'Heavy' seeding scenario 

The origin of the first massive black holes is still a subject 
of intense deb ate. Currently , there exist two f avored seeding 
mechanisms (|Haimanl |20ld : IVolonteril |2010h : either black 
hole seeds could form out of the remnants of massive Pop III 
stars (e.g. iMadau fe Reedl200ll ; iHeger fe Wooslevll2002l ) or 



during the direct core-collapse o f a low - angular momentum 
gas cl oud (e.g. 



2004 IVolonteri et 
Bellovarv et alj 120111 



it al.1 



Loeb fc Rasiol 1 1994 : Koushiappas et al] 
l200Sl : IVolonteri fc Starkl 1201 
In the first case the seeds are 
expected to have masses of M sco d ~ 100 — 600Mq ('light' 
seeding), while in the latter case, more massive seeds 
between M aee d ~ 10 5 — 10 6 Mq ('heavy' seeding) would be 
expected. The detailed physical processes, in particular of 
the direct core-collapse, are largely unknown. Unfortunately, 
observational constraints in the high-redshift universe are 
too weak to favor one of these models. Ho wever, future 
obser v ations of gravitational waves (LISA, 1 esana et al.l 
120051 : iKoushiappas fc Zentner 20061) o r planned X-ray 
missions (WFXT: ISivakoff et al] l2010l : iGilli et al] I2OI0I ; 
IXO), may have the ability to detect accreting black holes 
at z > 6, and thus, will be able to test these models of the 
first black holes. Moreover, due to the exponential growth 
of the black holes during accretion, it is also very difficult 
to use the local population of massive black holes to recover 
information about their original masses before the onset 
of accretion. For ins tance, in the t h eoret ical studies of 
IVolonteri et~ai1 (|2008h . and IVolonteril <|2010h . the different 
seeding mechanisms are investigated by following the mass 
assembly using Monte-Carlo merger trees to the present 
time. They find that both models can fit observational 
constraints at z — (e.g. the black hole mass-velocity 
dispersion relation or the black hole mass function), when 
light seeds form already at very early times (z = 20), or 
when heavy see ds evolve later on (z = 5 — 10). Furthermore 
in a study of iTanaka fc Haimanl (|2009h they use dark 
matter halo merger trees, coupled with a prescription for 
the halo occupation fraction and they show that ~ IOOA^q 
seed BHs can grow into 1O 6 M0 BHs at z ~ 6 without 
super-Eddington accretion, but only if they form in miniha- 
los at z ^ 30. In our baseline model, seed black hole masses 
of IOOMq were assumed, however, due to our adopted mass 
resolution, we obtain very little seeding before z ~ 10. 
Therefore, we also explore a heavy seeding scenario with 
M. iSee d = 10 5 Mq. The different seeding mechanisms will 
not affect the z = black hole mass and AGN population, 
as initial seed masses are compensated by gas accretion 
growth processes by orders of magnitude. Only the black 
hole distribution and QSO luminosity function at high 
redshifts will be influenced by this modification. Further- 
more, some studies have suggested that low mass dark 
matter halos may not be able to prod uce massive seeds 
lIMenci et al. 2008; IVolonteri et al] l201ll ). as the potential 
well might be too weak for collapse to occur. Therefore 
in the 'heavy seeding' model we additionally adopt a halo 
mass limit of 2 x 1O 11 M0, below which no black hole seeds 
are inserted. 

Furthermore 
2004 iPeng et al] 



many observations (e . g. IWalter et al 
20061: iMcLure et"al]|2006l:1Schramm et al 



2008) suggest that the black hole-to-bulge mass ratio was 



larger at higher redshifts than expected from the local black 
hole-bulge mass relation. This eventually implies that black 
holes were accreting more gas and thus, growing faster than 
the corresponding bulges at high redshifts than at lower 
ones. Therefore, besides assuming an evolving black hole- 
bulge mass relation (see the z-dependent F parameter in 
eq. [TJ we additionally adopt a larger scatter oy accr for 
the accreted mass onto the black holes at high redshifts. 



Origin of the anti-hierarchical growth of black holes 7 



This means that, when calculating the final black hole mass 
M. , final, a larger Gaussian distributed scatter is applied with 
a value of <7., a ccr = 0.6 for z > 4. At redshift z < 4 the orig- 
inal scatter value cr. iaccr = 0.3 is applied. We find that we 
still recover a tight relationship between BH mass and bulge 
mass at z — 0, in agreemen t with observation s . This is in 
agreement with the results of lHirschmann et all (|2010l ). who 
showed that a large scatter in black hole mass at fixed bulge 
mass (<r = 0.6 dex) at high redshift (z = 3) will decrease 
towards the observed present-day value due to mergers. 



3.5 Summary of Model Variants 

In the course of this study, we investigate the effects of 
the outlined modifications on the AGN/black hole evolu- 
tion one-by-one and in various combinations. We consider 
the following six different models: 

1. FID: Fiducial, standard accretion model (|3 . 1 f) 

2. VE: Varying sub-Eddington limit for the maximum ac- 
cretion rate / e dd H3.II & 13. 2[) 

3. DI: Additional accretion due to Disk Instabilities (|3.1I 
fe l33|> 

4. SH: Heavy Seeding mechanism with a Halo mass limit 

5. DISH: Disk Instabilities & Heavy Seeding mechanism 
with a Halo mass limit (|3. 1 13.31 &: [3T4 |) 

6. VEDISH: Best-fit model including a Varying sub- 
Eddington limit, Disk Instabilities & a heavy Seeding 
mechanism with a Halo mass limit (|3.11 13.21 13.31 & I3.4[> 



4 PROPERTIES OF NEARBY GALAXIES AND 
BLACK HOLES 

In Fig. [T] we compare different galaxy and black hole prop- 
erties from the FID, the DISH and the VEDISH model to 
observations of the local Universe. We do not explicitly show 
the predictions of the VE, the DI and the SH model sepa- 
rately as they do not result in a stronger deviation from the 
FID model than the DISH or the VEDISH model. 

The upper left panel in Fig. Q] shows the mo deled stel- 
lar ma ss fu nctions compa r ed to observations from lBell et al.l 
(|2003h and iPanter et ail <|2007h . The modifications in the 
DISH and the VEDISH model hardly show any variation 
from the FID model and thus, for stellar masses larger 
than 10 9 M© we obtain a reasonably good match to the ob- 
servational data. However, low- mass galaxies (< 10 9 Mq) 
are slightly over-predicted, a commo n feature of most 
(all) current semi-analyt ic models (e.g. I Bower et al.l I2006I : 
|Pe Lucia fc Blaizotll2007l) and still a subject of on- going re- 
search (e.g.lGuo et al.ll201ll : iBower et af]|201ll : IWang et al] 
l201ll : iMenci et all I2012T ). Stronger supernova feedback, a 
modified star formation law or a different cosmological 
model are currently considered as possible solutions for this 
problem. 

The present-day black hole mass function is depicted in 
the upper right panel of Fig. [1] Our model predictions are 
compa r ed with observationa l esti mates from Shankar et al.l 
l|2004f ). lMarconi et al.l ( |2004l ) and lShankar et all (|2009l 'l (see 
the review of IShankarll2009l for more details). We find rea- 
sonably good agreement of the SAM predictions with the 



observations for the whole black hole mass range. Deviations 
between the different SAM models are negligible, indicating 
that neither the growth of black holes by disk instabilities 
nor the limited accretion rates at low redshifts influence the 
black hole mass function significantly. In contrast to our re- 
sult, in many SAM stud ies an excess of very massive black 

iMalbon et all 



Fontanot et al. 200 



holes can be seen (e . g. ^ 

120071 : iMarulli et all I2OO8I : iFanidakis et al.l bold ), which 



might be ca used by too-effi c ient R adio-mode accretion as 
discussed bv lFontanot et all |201ll ). 

The present-day relation between black hole and bulge 
mass is shown in the lower left panel of Fig. [T] with our 
model pre dictions compared wit h the observational relation 
derived bv lHaring fc Rixl |2004l ). All of the model variants 
reproduce a tight relationship between black hole mass and 
bulge mass, which results from the self-regulated BH growth 
assumed in our model. The slight upturn at low bulge masses 
in the DISH and the VEDISH models is due to the heavy 
seeding mechanism, where by construction no black hole 
masses below 10 5 M Q can exist. 

Finally, the lower right panel of Fig. [1] shows the ac- 
tive black hole mass function predi cted by our model, com- 
pared with observat i onal d ata from I Greene et all (|2010l ) and 
ISchulze fc Wisotzkil (|2010l ). Note that — consistent with 
these observational studies — we define "active" BH here 
as having a bolometric luminosity greater than 10 43 ' 5 erg/s. 
For M. < 1O 8,3 M0, the accretion due to disk instabili- 
ties (DISH/ VEDISH models) increases the fraction of active 
black holes by almost one order of magnitude compared to 
the FID model. In contrast, for M. > 10 8 ' 3 M© the limited 
gas accretion in the VEDISH model reduces the active frac- 
tion of massive black holes compared to the DISH and the 
FID model. Overall, the VEDISH model provides a better 
match to the observational data than the FID or the DISH 
models. 



5 GALAXY AND BLACK HOLE PROPERTIES 
AT HIGHER REDSHIFT 

Fig. [2] illustrates the stellar mass functions at different red- 
shifts (z = 0.5, 1, 2, 3, 4, 5) as indicated in the legend of each 
panel. We compare our SAM predict ions with the results 
of d i fferent observational studies 
|200S|; iBundv et al l 120051; iDrorv et 



s |pj 

lame 



Perez-Gonzalez et al 
2004; iFontana et al 



120061 : iMarchesini et al.ll2007l : lllbert et al.ll2010l ). At all red- 
shifts, there is no significant difference between the dif- 
ferent model variants. However, compared to observations, 
the high mass end is under-predicted, while the low- 
mass end is over-predicted by the SAMs. This discrep- 
ancy is again a well-known probl em ( Fontanot et al.l [20091 ; 
IMarchesini fc van Dokkuml 1200 il ). iFontanot et alj (|2009l ) 
sho wed this for t he M organa, Munich and the S08 model, 
and iGuo et all (|201ll ) found the same for the latest ver- 
sion of the Munich model. The discrepancy at the high- 
mass end may be related to systematic errors or scatter in 
the photometric s t ellar m ass estimates from observations 
— IFontanot et al. I (HH) showed that when stellar mass 
errors of 0.25 dex were convolved with the model results, 
the SAMs agreed reasonably well with the available stellar 
mass function c o mpila tions at least to z ~ 3. Moreover, 
ISomerville et al.l (|201ll ) showed that the SAM presented 
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Figure 1. Comparison of SAM results to observations for black hole and galaxy properties of the local Universe. The solid lines are the 
FID (green) and the VEDISH (red) model, the dashed, orange lines th e DISH mode l. Upp er lef t panel: the stella r mass function. The 
black solid line and the black symbols show observational results from iPa nter et al.l <l2007f) a nd Bell ct al. (2003), respectively. Upper 
right panel: the black hole mass function. Observational estimates dMarcom et al. 1 12004 IShankar et al. 1 12004 2009) are illustrated by the 
black solid lines, the grey shaded areas and the black symbols. Botto m left panel: black hole-bulge mass relation. The black solid line and 
symbols correspond to the observed relation bv lHaring fc Rixl <2004) and the grey shaded area show the l-cr scatter of the observational 
data. Bottom right pan el: the mass function of ac tive black holes with bolom etric luminosities above Lbol < 10 43 5 erg/s. Open symbols 
show observations from [Greene et al.l ll2010h and lSchulze &; Wisotzk Disk instabilities increase the number of active low mass 

black holes and the gas-dependent Eddington-limit decreases the number of active high mass black holes. 



here agrees with the observed rest-frame A'-band luminos- 
ity function at the bright end up to z ~ 3. The excess of 
low mass galaxies at high redshift is also seen in numerica l 
cosmological hydrodynamic simulations (iDave et al.ll201ll ). 
showing that this problem is not peculiar to SAMs. Instead, 
it may be an indication that the star formation or super- 
nova feedback recipes that are commonly adopted in both 
SAMs and numerical simulations require revision (Caviglia 
& Somerville, in prep.), or that the underlying cosmological 
model differs from the Cold Dark M atter paradigm (e.g., 
warm dark matter: [Menci et al.ll20l3 ). 

Fig. [3] shows the black hole mass function at different 
redshift steps (z = 0,0.5,1,2,3,4,5,6) for the FID, DISH, 
and VEDISH models. We show the observational estimates 
at z — to guide the eye. At redshifts z < 3, the dif- 
ferent model assumptions in the DISH and the VEDISH 
models do not influence the evolution of the black holes. 



However, turning to higher redshift z > 3, the main differ- 
ence between the FID model and the DISH/ VEDISH mod- 
els is the larger number density of black holes more massive 
than M, > 1O 6 M0. This can be explained by the 'heavy' 
seeding scenario and the large scatter in the accreted gas 
mass onto the black hole at z > 4. This leads to larger 
black hole masses and faster growth at early redshift than 
in the FID model. Towards lower redshift, however, this ef- 
fect dissappears as the subsequent growth by gas accretion 
overcomes the seed black hole masses by orders of mag- 
ni tude. This trend is eve n more pronounced in the model 
of iFanidakis et ail (2010), as they allow the black holes to 
accrete at super-Eddington rates: e.g. at redshift z — 6, 
black holes with M. = 1O 6 M0 have a number density of 
log$ = —2.7 Mpc -3 dex" 1 , whereas we obtain a number 
density of only log$ = -4 Mpc" 3 dex -1 in the VEDISH 
model. 
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Figure 2. Redshift evolution of the stellar mass function (2 = 0.5,1,2,3,4,5). The solid lines show the FI D (green) and the VEDISH 
(red) model and the orange dashed lines the DISH model. The black da shed lines depict obse r vation s from|l lbert et al.l 1 2010TI and the 
symbo l s correspond to obser v ations from a set of obse rvational studies ( Perez-Gonzale z et al. I 120081 ); iBundv" et al. I <200a) : lDrorv" et aU 
[|2004h : iFontana etaD <|2006l l: iMarchesini et alj j2007h ). As in most SAMs, the number of low mass galaxies is overestimated (at all 
redshifts), while the number of high mass galaxies is underestimated at high z. 



The evolution of the active fraction of black holes is 
shown in Fig. [4] The three panels correspond to the dif- 
ferent models (FID: upper panel, DISH: middel panel and 
VEDISH: lower panel), where colored lines illustrate the 
model results at different redshifts (z — 0, 0.5, 1, 2, 3, 4, 5, 6). 
We consider only AGN with bolometric luminosity larger 
than 10 43 ' 5 erg/s. The grey symbols show the local obser- 
vations. We find that the evolution of the active black hole 
fraction varies from model to model, implying that our mod- 
ifications to the recipes for BH formation and evolution have 
a significant influence on the active fraction at all redshifts. 
Comparing the DISH model with the FID model we can see 
two effects: At high redshifts 2 3, the number of active 
black holes with masses between 10 6 < M. < 1O 8 M0 is 
greatly increased due to the heavy seeding mechanism and 
the large scatter in the accreted gas mass. At low redshifts 
2^1, the number of active black holes with masses be- 
low 10 8 ' 3 Mq rises as a consequence of the additional gas 
accretion due to disk instabilities. Furthermore, as already 
seen in Fig. [1] the limited accretion rate in the VEDISH 
model reduces the fraction of AGN with massive black holes 
(> 10 8 ' 3 M o ) at z < 0.5. 



6 NUMBER DENSITY EVOLUTION OF AGN 

The different panels in Fig. [S] show the redshift evolution of 
the AGN number densities as a function of the bolometric 
luminosity, for the six different models. In this section, our 



SAM pre dictions are compared to the observational compila- 
tion from I Hopkins et al.l (|2007l ). In their study, they convert 
the AGN luminosities from different observational data sets 
and thus, from different wavebands (emission lines, NIR, 
optical, soft and hard X-ray) into bolometric ones. They 
assume a luminosity dependence of the obscured fraction 
(the less luminous the more obscured) and the same num- 
ber of Compton-thick (Nh > 10 24 cm -2 ) and Compton-thin 
(10 23 cnT 2 < N H < 10 24 cm" 2 ) AGN. However, there are 
many aspects of the obscuration corrections that are still 
being vigorously debated. Some recent studies suggest that 
the obsc ured fraction is dependent on bo th luminosity and 
redshift l|Hasingej|200i : iFiore et al.ll2012l). in contrast with 
the non-redshift dependent model of I Hopkins et al] (2007). 
There are also uncertaintie s surrounding t he du st correc- 
tion for the UV luminosity; Hop kins et al.l (2007) compute 
the amount of dust (and therefore extinction), by adopt- 
ing an Nh distribution from X-ray observations, and a 
Galactic dust-to-gas ratio. However, it has been shown that 
AGN absorbers do not hav e a Galactic dust to gas ratio 
|Maiolino et al.ll200ll . 12004 ) ■ The result is that they prob- 
ably over-estimate the extinction, which might result in 
slightly higher luminosities for the optically selected quasars 
(F. Fiore, personal communication). Because of these uncer- 
tainties, we both com pare the obscuration- corrected obser- 
vational compilation of lHopkins et al.l (|2007l ) with our unob- 
scured model predictions, and in Section \7. 21 we attempt to 
correct our model predictions for obscuration and compare 
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Figure 5. Number densities of AGN versus redshift for the six different SAM models (FID, VE, DI, SH, DISH and VEDISH) as indicated 
in the legend. Different colors illustrate different bolomctric luminosity bins: red: 47.5 < log(Lt, i), yellow: 46.5 < log(Lt, i) < 47.5, green: 
45.5 < log(Lboi) < 46.5, light blue: 44.5 < log(L b ol) < 45.5, dark blue: 43.5 < log(L b ol) < 44.5, black: 42.5 < log(L bo i) < 43.5. Solid 
lines and open squares show t he corresponding mod el predictions; dashed lines and stars together with grey shaded areas indicate the 
observational compilation from Hop kins et al.l 1(2007). While the FID model shows the opposite of "downsizing" behavior, with luminous 
objects forming late and low-luminosity objects forming early, we obtain a fairly good match to the observations for the VEDISH model. 



Origin of the anti-hierarchical growth of black holes 11 




6 7 8 9 10 

log(M./M ) 



Figure 3. Redshift evolution of the black hole mass function. 
The colored lines illustrate the SAM results of the FID (upper 
panel), DISH (middle panel) and the VEDISH (low panel) model 
at different redshifts. T he grey lines and symbols correspond to 
obser vational estimates l|Marconi et al J2004l ;l Shank ar et al .120041 . 
2009) at z = 0. Massive black hole seeds significantly increase the 
number of black holes with masses between 10 6 < M. < 10 s Mq 
at z > 4. 



with recent soft and hard X-ray measurements of the AGN 
luminosity function. 

The upper left panel in Fig. [S]shows the FID model. The 
number densities at the peak of each luminosity bin are in 
quantitative agreement with the observations implying that 
the FID model reproduces the correct order of magnitude 
of AGN number densities in the different luminosity bins. 
However, the observed time evolution of the peaks of the 
different luminosity classes are not correctly predicted by the 
FID model, which shows the typical "hierarchical" behavior 
in which the number of low-luminosity objects peaks earlier 




log(M./M ) 

Figure 4. Redshift evolution of the black hole mass function 
of active black holes with bolometric luminosities larger than 
^bol > 10 43 5 erg/s. Colored lines illustrate the SAM results of 
the FID (upper panel), the DISH (middle panel) and the VEDISH 
(lower panel) model at different redshifts. For comparison, the 
grey symbols illus trate observations of the active black hole mas s 
function at z = ijGreene et al.ll20ld : ISchulze fc WisotzkilEuTol ) . 
Massive seeds increase the number of active black holes at high 
redshift and disk instabilities increase the number of active low 
mass black holes at low redshift. 



than the number density of higher luminosity objects. We 
can summarize this problem more quantitatively as follows: 

• z < 2: over-prediction of AGN with log(Lboi) > 46 

• 2 < 2: under-prediction of AGN with log(Lboi) < 46 

• z > 3: under-prediction of AGN with log(Lboi) > 46 

• z > 3: over-prediction of AGN with log(Lboi) < 45 

As the fiducial model reproduces the black hole mass 
function at z=0, we can assume that there are the correct 
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Figure 6. Bolometric quasar luminosity functions at different redshifts. Black stars show the observational compilation of lHopkins et all 
(2007). The green solid lines correspond to the output of the FID model, the purple, dashed lines show the VE model, the light blue, 
dashed lines the DI model and the dark blue, dashed lines the SH model. 



number of black holes. Therefore the first point suggests that 
at low redshifts either too high a fraction of massive black 
hole are accreting or these massive black holes are accreting 
at rates that are too high. Moreover, assuming that activity 
is triggered by merger events implies that the natural de- 
crease in the major merger rate is not sufficient to produce 
the observed steep decline in the AGN number densities. 
It has been shown that the galaxy mer ger rate predicted 
by the S08 mod els matches observations (|Jogee et al.ll200ci 
lLotz et afl|201lf ) so this is not likely to be the cause of the 
discrepancy. The low number densities of moderately lumi- 
nous AGN may indicate, however, that the AGN activity 
might not only be triggered by merger events, but also by 
secular evolution processes. The deviations at high redshift 
may be a consequence of massive black holes forming too late 



in the SAM as well as possibly the non-redshift-dependent 
dust obscuration correction. The excess of moderate and 
low luminosity AGN at high redshift may be partly due to 
the similar over-prediction of low-mass galaxies in the SAM, 
which we already discussed in Section [5] 

We experimented with whether we could achieve bet- 
ter reproduction of the observed downsizing of BH activ- 
ity by just modifying the values of some free parameters in 
the FID SAM, without changing any physical ingredients. 
For example, we tested the effect of varying the strength 
of Supernova and r adio-mode feedback (as suggested in 
iFontanot et al.|[2006l ). We find that for stronger supernova 
feedback (doubling the normalization of the mass loading 
factor for SN-driven winds £sn) we achieve a decrease in the 
number density of AGN at all redshifts, resulting in over- 
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Figure 7. The same as in Fig. [6] but for different SAM models: the green solid lines correspond to the output of the FID model, the 
orange, dashed lines show the DISH model and the red, solid lines illustrate the results from the VEDISH model. For the VEDISH model, 
a reasonably good agreement with observations is obtained. 
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all a worse match to the observations. Even if the number 
density of moderately luminous objects (Lboi < 10 45 erg/s) 
at high redshifts decreases (towards the observational data), 
this is not sufficient to achieve a reasonably good match to 
observations in this range. Increasing the strength of the 
radio-mode feedback reduces the number of luminous AGN 
(Lboi > 10 45 erg/s) between redshifts < z < 4, again 
resulting in worse agreement with the observational com- 
pilation than for the FID model itself. If we assume 'halo 
quenching' instead of the radio-mode feedback model, i.e. 
no cooling is allowed above a certain threshold halo mass 

o,thres — 

1O 12 M ), we obtain a decrease of the number 
of luminous AGN, but again at all redshifts, and not redshift 
dependent as observed. To increase the number density of 
moderately luminous objects at low redshifts, we varied the 



timescale (tq) in th e power-law decline phase of a quasar 
episode. A study bv lMarulli et al.l (|2008l ) has shown that a 
power-law decline growth phase in their quasar mode does 
increase the number density of moderately luminous AGN at 
low redshift, resulting in a better agreement with the obser- 
vations. However, within our study it is found that varying 
t<3 is not sufficient to match the observations. Therefore, we 
can conclude that downsizing cannot be reproduced solely 
by varying the free parameter values of our FID model. In- 
stead, we now present the influence of the additional modi- 
fications for black hole formation and growth as outlined in 
Section [3] 

The result for the VE model (i.e. assuming a sub- 
Eddington accretion limit dependent on the cold gas frac- 
tion) is illustrated in the upper right panel of Fig. [5] For 
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bolometric luminosities larger than 10 46 erg/s, the cold gas 
fraction dependent sub-Eddington limit clearly reduces the 
number density of AGN. We find that we are able to re- 
produce the observed steep decline in the number densities 
of bright AGN at z < 2 reasonably well in this way. This 
supports the idea that the cold gas content of a galaxy may 
regulate the efficiency of black hole accretion, in particu- 
lar the maximum accretion rate that can be reached in a 
merger. We may speculate that low cold gas densities lead 
to smaller viscosities so that it takes longer for the gas to 
lose its angular momentum and thus, to be accreted onto 
the black hole. Our results imply that even when massive 
black holes experience major mergers, if the gas fractions in 
their host galaxies are low they may not produce luminous 
AGN because they never accrete at close to the Eddington 
rate. 

The middle left panel of Fig. [5] shows the results of the 
DI model, assuming an additional BH accretion mode due to 
disk instabilities. For a luminosity range of 43 < log(Lbol) < 
45, the number of AGN is increased, resulting in a rea- 
sonably good match with the observational compilation for 
2 < 1.5. For the lowest luminosity bin, however, the number 
of AGN is now over-predicted. Nevertheless, this additional 
accretion mechanism seems to play an important role for 
triggering the activi ty of faint AGN (co nsistent with obser- 
vational results e.g. Sa lucci et al.l feoOO'). In strong contrast 
to our model, black hole ac cretion due to disk in stabilities 
in the Galform- model (e.g. lFanidakis et al.ll2010r i provides 
the major contribution to AGN number densities for all lu- 
minosities a nd at all redshifts. I n the MuNiCH-models a s 
presented bv lMarulli et al.l (120081 ) and iBonoli et all l|2009l ). 
black hole accretion due to disk instabilities is not accounted 
for at all, but they still slightly under-predict the number 
of moderately luminous AGN, even in their best-fit model. 
This further supports the need for BH accretion driven by 
secular evolution processes, in addition to mergers. 

The effect of a heavy seeding mechanism together with 
a halo mass limit for black hole formation (SH model) is il- 
lustrated in the middle right panel of Fig. [5] The number of 
bright AGN at high redshift is increased and can match the 
observational data. This is a consequence of large seed black 
hole masses and a large scatter in the accreted gas mass. As 
we cap BH growth at the Eddington rate, having more mas- 
sive seeds means that these early black holes can grow faster, 
leading to a larger number of massive and active black holes 
at early times than in the FID, VE and DI models. Our re- 
sult indicates that black holes probably have to undergo a 
phase of very rapid growth at high redshifts (z > 5), even 
if it is still unknown whether and how such massive seed 
black holes can form out of direct core-collapse or whether 
less massive seeds have to accrete at super-Eddington rates. 
However, we find that assuming even more massive black 
hole seeds with masses of A/. iSCC d = 1O 6 M0 results too few 
moderately luminous AGN at high redshift (z ~ 5) in our 
model. The halo mass limit for seed black hole formation also 
reduces the number density of faint AGN, as black holes in 
low-mass halos are not allowed to form and to accrete gas. 
However, even if this second effect results in better agree- 
ment with the observational data, it does not seem to be 
fully sufficient for reproducing them. This might indicate 
that the dependence of seed mass on halo mass is more com- 
plex than we have assumed in this simple model, or might 



be be due to redshift-dependent obscuration, which has not 
been accounted for in the observational comparison that we 
are using here (see Section [7.211 . 

Finally, the combination of the individual modifications 
which have been discussed so far is presented by the DISH 
and the VEDISH model (lower left and lower right panel of 
Fig. [5]). We find that the changes in the AGN number densi- 
ties seen for the individual modifications sum in a straight- 
forward way, without significantly influencing each other. 
Thus, the VEDISH (= "best-fit" ) model is able to reproduce 
the observed downsizing trend fairly well, and can predict 
the correct time-evolution of the peaks of the luminosity de- 
pendent number density curves. However, even in our best- 
fit model, the number of faint AGN is still over-predicted 
at redshifts between 4 > z > 2. This might be due to the 
difficulty of detecting these objects at high redshifts, either 
because they may not be easily recognizable as AGN, or 
because of obscuration. 



7 THE AGN LUMINOSITY FUNCTION 
7.1 Bolometric luminosities 

The effects of our individual modifications for black hole 
growth and the final success of the VEDISH model can also 
be explicitly seen in Figs. [6] and [7] where the bolometric 
AGN luminosity function (AGN LF) is plotted at different 
redshifts, and again compared with the observational com- 
pilation from [Hopkins et al.l (2007). Fig. [6] shows that the 
DI model raises the low-luminosity end by about one order 
of magnitude in AGN number densities, while the VE mod- 
els lowers the high-luminosity end by more than two orders 
of magnitude at z ^ 1.5. The SH model changes the AGN 
number densities at z ^ 3 by lowering the low-luminosity 
end and raising the high-luminosity end, each of them by 
about one order of magnitude. The cumulative effect of the 
separate alterations is shown by the DISH and the VEDISH 
model (Fig. [7|). While the DISH model still fails to repro- 
duce the high luminosity end at low redshift, the VEDISH 
model represents a fairly good match with the observational 
data for the whole redshift range. 



7.2 X-ray Luminosity Functions 

We now compare our model predictions with recent obser- 
vational determin ations of the AGN LF in the hard and 
soft X-ray bands jHasinger et al l 120051 ; lEbrero et all 120091 ; 
lAird et al.ll20ld ; iFiore et al.ll2012T ). In contrast to the previ- 
ous section, we do not attempt to correct the observations 
for obscuration, but instead apply an obscuration correction 
to our models. We convert the modeled, bolometric lumi- 
nosities into hard and soft X-ray luminosities (0.5 — 2 keV 
and > 2 keV) using the bolometric correction according to 
iMarconi et al.l (|2004T ). In their study, the hard and soft X-ray 
luminosities Lhxr, Lsxr are approximated by the following 
third-degree polynomial fits: 



log(i H xR/iboi) = -1.54 - 0.24£ - 0.012£ 2 + 0.0015£ 3 (14) 
log(L S xR/iboi) = -1.65 - 0.22£ - 0.012£ 2 + 0.0015£ 3 (15) 

with £ = log^boi/^©) — 12. These corrections are derived 
from template spectra, which are truncated at A > 1 fim in 
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Figu re 8. Soft X-ray AG N luminosity f u nction for different redshift ranges. Black stars and circles show the observational data from 
lHasinger et alj HqoH) and lEbrero et all (120091 1. respectively. The solid green, dashed orange and solid red lines correspond to the 
results of the FID, the DISH and the VED I SH m odel and the blue solid lines illustrate the results of the VEDISH model with an 
obscuration correction according to Hasinger (2008). The VEDISH model including the obscuration correction is in fair agreement with 
the observational data. 



order to remove the IR bump and which are assumed to be 
independent of redshift (therefore the resulting bolometric 
corrections are also assumed to be redshift independent). 
Additionally, we apply a correction for obscuration to the 
model l uminosities, as suggested by several observational 
studi es l|Ueda et ail 120031: iHasingeif |2004| ; lEa Franca et al.l 
2005), in which it has been shown that the fraction of ob- 
scured AGN is luminosity dependent and decrea ses with in- 
creasi ng lu minosity. Whi l e olde r studies such as lUeda et al.l 
(|2003h and ISteffen et all (|200St ) did not find a clear depen- 
dence of obscuration on redshift , several recent observational 
studi es (|Ballantvne et al.l 120061 ; iGilli et all 120071 ; lHasingerl 
2008) propose a strong evolution of the obscured AGN pop- 
ulation (with the relative fraction of obscured AGN increas- 
ing with increa sing redshift). Here, we follow the study of 
lHasingerl |2008), where they compare the same AGN in both 
the soft and hard X-ray band so that they can derive an ap- 
proximation for the obscured fraction in the soft X-ray band. 
The obscured fraction is then given by this equation: 



fa 



-0.281(log(L LX a) - 43.5) + 0.279(1 + z) a . 



(16) 



where they find that a value of a — 0.62 provides the best 
fit to their observational data. By calculating the obscured 
fraction of AGN in the soft X-ray band, we can model the 
visible fraction of AGN / v i s = 1 — /obsc and thus, the visible 
number density of AGN in the soft X-ray range is given by: 

$vis(£sXB.) = /vis X $total(isXR.) (17) 

Fig. [8] shows the soft X-ray luminosity function for 



different redshift ranges predicted by our FID, DISH 
and VEDISH model s, co mpared with o bserv ations from 
lHasinger et"afl |2005l ) and lEbrero et all (|2009l ). We show 
the model predictions with and without the obscuration 
correction described above. For the high-luminosity end 
(isxR > 10 45 erg/s), obscuration does not influence our 
results and thus, one can see the same trends as already dis- 
cussed in Fig. Turning to the low-luminosity end at low 
redshifts (z ^ 0.8), the FID model matches the observational 
data, while the DISH and VEDISH model overproduce the 
number density of AGN. However, considering obscuration 
effects leads to better agreement of the VEDISH model with 
the observed number densities than the FID model. At high 
redshift (z > 0.8), both the FID and the VEDISH model 
over-predict the number of moderately luminous AGN by 
almost the same amount, and both models including obscu- 
ration can achieve a fairly good agreement with the obser- 
vations. Overall, the VEDISH model including AGN obscu- 
ration can predict the observed soft X-ray luminosity func- 
tion reasonably well, support ing the adoption of a redshift- 
dependent obscured fraction. iFanidakis et al.l l|2010t ) use the 
same bolometric conversion for calculating soft X-ray lu- 
minosities and the same dust obscuration correction as in 
our study. They show that they are able to match the soft 
X-ray luminosi t y fun ctions from the observational study of 
lHasinger et~ai1 J2005h as well. 

Fig. [9] illustrates the hard X-ray luminosity func- 
tions predicted by our models for diffe rent redshift ranges , 
compared with observational data (Eb rero et al.l 120091 . 
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Figur e 9. Hard X-ray AGN luminosity functi on for different redshift ranges. Black symbols show the observational data of lEbrero et al.l 
j200Sh; lAird et aU <201CD ; iFiore et all <2012l) and data from the SDSS (optical luminosity is converted into X-ray luminosity as in 
IFiore et alTl2012l) '). The solid green, dashed orange and solid red lines illustrate the output of the FID, the DI SH and the VE DISH 
model. The blue solid lines illustrate the results of the VEDISH model with an obscuration correction according to Hasinecr (2008). The 
latter is able to predict the observed hard X-ray LF reasonably well. 
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lAird et all l20ld , iFiore et all |2012| . iFiore et all |2012| ). The 
VEDISH model is able to reproduce the high-luminosity end 
pretty well, whereas it over-predicts the number of AGN 
at the low-luminosity end at all redshifts (with a larger 
discrepancy at higher redshift). One possible explanation 
might be again due to obscuration as even 2 — 10 keV X- 
ray surveys might miss a significant fraction of moderately 
obscured AGN ( 25% at N H = 10 23 c m" 2 ) and nearly all 
Compton-thick AGN ( N H > 10 24 cm~ 2 . lTreister et al.ll2004l : 
Balla ntvne et a l. 2006). From fits to the cosmic X-ray back- 
ground. lGilli et all (|2007h predict that both moderately ob- 
scured and Compton-thick AGN are as numerous as un- 
obscured AGN at luminosities higher than log(Lo.5-2k c v) > 
43.5[ergs/s], and four times as numerous as unobscured AGN 
at lower luminosities (log(Lo.5-2k c v) < 43.5[ergs/s]). For 
this reason, we made the very simplified assumption that 
the fraction of obscured AGN in the hard X-ray band is 
the same as in the soft X-ray band. This results in a fairly 
good match to the observed low-luminosity end at all red- 
shifts, except around z ~ 1 where the obscuration-corrected 
model under-predicts the number of faint AGN. Therefore, 
our model results are consistent with the existence of an ob- 
scured fraction of AGN in the hard X-ray band (Compton- 
thick AGN) that is of the same order of magnitude as 
the obscured fraction in the soft X-ray band. Interestingly, 
iFanidakis et al. (2010) show that they are able to match 
the hard X-ray luminosity functions from an observational 
study of lUeda et aTJ l|2003l ) without including any obscura- 
tion effects. However, their predictions are only illustrated 
for a comparatively small redshift range of 0.2 < z < 1.6 and 
for hard X-ray luminosities larger than Lhxr > 10 42 erg/s, 
where our unobscured VEDISH model predictions are also 
in good agreement with the observational data. 



8 EDDINGTON RATIO DISTRIBUTIONS 

Having assessed our models by comparing with the ob- 
served AGN LF, we now examine the consequences of our 
modifications of the black hole growth prescriptions on the 
Eddington ratio distribution and on the AGN luminosity- 
black hole mass plane. In Fig. [TU] we show the redshift 
evolution of the Eddington ratio distributions for the FID 
(upper panel), the DISH (middle panel) and the VEDISH 
model (lower panel). Different colors indicate different red- 
shift steps (z — 0,0.5,1,2,3,4,5,6). Consistent with our 
previous investigation, and typical limits of observational 
samples, we consider only AGN with bolometric luminosi- 
ties larger than 10 43 erg/s. In all models, the fraction of 
AGN that are are accreting at smaller Eddington ratios in- 
creases strongly with decreasing redshift. This is because at 
later times, black holes spend less of their time in the first 
regime, i.e. accreting at the Eddington rate, but mainly re- 
side in the second regime, the power-law decline dominated 
"blowout" accretion phase. These black holes are relics from 
an earlier, more active phase with higher accretion rates. 
While, however, in the FID model the Eddington ratios are 
still peaking between / c dd =0.1 — 1 at all redshifts, in the 
DISH and the VEDISH models the peaks of the distribu- 
tion curves are clearly shifted towards smaller Eddington 
ratios with decreasing redshift. This is mainly caused by 
the additional black hole accretion due to secular evolu- 
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Figure 10. Eddington ratio distributions for the FID, DISH 
and the VEDISH models (upper, middle and lower panel panel, 
respectively). Different colors correspond to different redshifts 
(z = 0.5,1,2,3,4,5,6). Note that here we show only accreting 
black holes with bolometric luminosities larger than 10 43 erg/s. 
The modifications in the VEDISH model cause the peaks of the 
distributions to shift towards smaller Eddington rati os with time 
in qualitative agreement w i th observational studies jV cstcrgaard 
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tion processes. At 2 = 6 the distributions of the DISH and 
VEDISH models peak at / e dd ~ 0.1, while at z = the 
peaks are located around / c dd ~ 0.01. Thus, the majority of 
AGN at z — are not radiating at or near the Eddington 
limit anymore, whi ch is in qualitative agreement with ob- 
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tion (using broad-line quasars between 0=1 — 4) peaks 
at an Eddington ratio of / c dd = 0.05. Furthermore, the 
large seed black hole masses and the large scatter in the 
accreted gas mass in the DISH and VEDISH models lead 
to a significantly larger number of accreting black holes 
with Eddington ratios between 0.01 < / c dd < 1 at z > 4 
than in the FID model. Finally, the VEDISH model addi- 
tionally shows that as a consequence of the sub-Eddington 
accretion rate limit (dependent on the cold gas fraction), 
the number of black holes accreting close to the Eddington 
limit at z = is reduced by about one order of magni- 
tude compared to the DISH model. Altogether, we find that 
the downsizing behavior seems to imply that the peaks of 
the Eddington ratio distributions are shifted towards smaller 
Eddington ratios with decreasing redshift, and the number 
of black holes accreting close to the Eddington-rate is low at 
2 = (\og(dN/df cdd ) = -6.5 Mpc -3 dex" 1 ), while at high 
redshifts a large number of black holes (log(dAT/df e dd) = 
—5.2 Mpc -3 dex -1 ) are accreting within a broad range of 
Eddington ratios 0.01 < / e dd < 1. 



9 LUMINOSITY-BLACK HOLE MASS-PLANE 

Fig. [TT] shows the bolometric AGN luminosity-black hole 
mass- plane at z — (left column), z — 2 (middle column) 
and 2 = 5 (right column) for different SAMs (first row: 
FID model, second row: VE model, third row: DISH model 
and fourth row: VEDISH model). The circles in each panel 
correspond to our model AGN with a bolometric luminosity 
cut of Lboi > 10 43 erg/s and they are color coded according 
to their bolometric luminosity as defined in Fig. [5] We also 
show lines indicating accretion at the Eddington rate, an 
Eddington ratio of / c dd =0.1 and an Eddington ratio of 
/cdd = 0.01. In addition, we sho w the observational l imit i n 
the Lboi — M.-plane according to lSteinhardt fc Elvisl (2010). 

At 2 = 5 (right column of Fig. [TT), in the FID and VE 
models, black holes within a mass range of 10 4 - 7 < M. < 
10 s Mq are active, where the luminosity is almost linearly 
correlated with black hole mass, indicating accretion at or 
very close to the Eddington rate. As a consequence of the 
heavy seeding mechnism and of the halo mass limit for black 
hole formation, the mass range of active black holes is shifted 
towards larger values 10 5 ' 5 < M. < 1O 9 ' 5 M in the DISH 
and the VEDISH models. As we showed previously, in or- 
der to reproduce the observed downsizing trend these bright 
AGN (L bo} > 10 46 - 5 erg/s) need to appear already at these 
early times. Turning to z = 2 (middle column of Fig. [TTj) . we 
see that in all models the number of AGN with massive BH 
(Mm > 1O 9 M ) has incre ased relative to z = 5. Co mpared 
to the accretion limit of ISteinhardt fc Elvisl l|201Ch . all of 
our model results are in acceptable agreement, although we 
predict some massive BH accreting at or very close to the 
Eddington rate which strictly speaking are not 'allowed' ac- 



cording to the lSteinhardt fc Elvisl |2010l ) results. Also for all 
models, the relation between black hole mass and bolomet- 
ric luminosity becomes much broader — black holes with 
masses Mm « 10 8 Mq can now also power moderately lu- 
minous AGN with Lboi ~ 10 43 erg/s as they are accreting 
with Eddington ratios below / c dd < 0.01. The probability for 
black holes with M, > 1O 7 M0 to accrete at Eddington ratios 
below /odd = 0.01 is even higher than to accrete at larger Ed- 
dington ratios. This is due to the power-law decline accretion 
phase that black holes are experiencing: massive black holes 
powering moderately luminous AGN are remnants of former, 
high-luminous AGN. In the DISH and VEDISH models, we 
can additionally see the effect of black hole accretion due to 
disk instabilities: the number of AGN with black hole masses 
of 1O 7 M < Mm < 10 s Mq and with luminosities between 
10 43 < Lboi < 10 45 erg/ s is significantly increased, accreting 
with Eddington ratios around / e dd ~ 0.01. As Seyfert galax- 
ies are mainly spiral galaxies with black hole masses in the 
range of lO 7 A/ < M. < 1O 8 M0 and are — compared to 
quasars — only moderately luminous, disk instabilities seem 
indeed to provide one of the most important trigger mecha- 
nisms for their nuclear activity. Moreover, in the DISH and 
VEDISH models, no black holes below 1O 5 ' 5 M are active, 
in contrast to the FID and the VE models. This is due to 
the halo mass limit for black hole formation. 



Finally, at redshift 2 = (left column of Fig. Ill[) . the 
number of actively accreting black holes is significantly re- 
duced compared to z = 2. The DISH and the VEDISH mod- 
els show again the effect of the additional accretion chan- 
nel due to disk instabilities, increasing the number of mod- 
erately luminous AGN, with typical black hole masses of 
Seyfert galaxies (1O 7 M < M, < 10 s Mq). While, however, 
in the FID and the DISH models, bright AGN can exist 
with massive black holes accreting at and close to the Ed- 
dington rate, in the VE and VEDISH models the limited 
accretion regulated by the cold gas fraction suppresses the 
appearance of the se bright AGN. Compared to the observed 
accretion limit of ISteinhardt fc Elvisl (J2010) , our results in 
the VE and VEDISH model are in very good agreement with 
the observed limit. This strongly suggests that the depen- 
dence of the accretion rates on the cold gas content might 
provide a possible physical origin for the observed accretion 
limit and thus, the reduced number of bright AGN at low 
redshifts. 



Interestingly, in the study of iFanidakis et all l|201Cf ). 
they obtain almost no evolution of their luminosity-black 
hole mass relation within a redshift range 0.5 < z < 2. Black 
holes above M, > 1O 9 M always accrete below / c dd = 0.01 
(ADAF regime), as in our VE and VEDISH models. How- 
ever, in clear contrast to this work, in this redshift range, 
they still have super-Eddington accretion, in particular for 
black hole masses between 1O 8 M < M. < 1O 9 M (re- 
sulting in Lboi > 10 46 erg/s). This seems to be, how- 
ever, in contrast to several observational studies, such as 



ISteinhardt fc Elvisl (|2010T ). which do not find black holes 
with masses M, > lO 7 Af accreting near or at the Edding- 
ton limit at low redshifts z = 0.2 — 0.4. 
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Figure 11. Bolometric luminosity versus black hole mass at z = (left column), 2 = 2 (middle column) and 2 = 5 (right column). 
Different rows correspond to the FID (first row), the VE (second row), the DISH (third row) and the VEDISH (fourth row) model. The 
open circles illustrate actively accreting black holes at each redshift. The different colors correspond to different bolometric luminosity 
bins as defined in Fig. [5] The red solid lines always illustrate accretion at the Eddington rate / e< jd = 1 , while the red dashed and dotted 
dashed lines depict Eddin gton ratios of / c dd = 0-1 an d fedd = 0.01, respectively. The green solid lines show the sub-Eddington limit from 
ISteinhardt fc ElviJ ll2010h . We find that the observed upper envelope in Eddington ratio is reproduced well by our adopted gas fraction 
dependent accretion limit. 



10 SUMMARY AND DISCUSSION 

In this study, we have used the semi-analytic model of S08 
to study the origin of the observed anti-hierarchical be- 
havior of BH activity within the framework of hierarchi- 
cal structure formation. In the original S08 SAM (FID), all 
AGN are merger-driven, with the prescription for BH ac- 
cretion and the light curve models based on hydrodynamic 
binary merger simulations. In these models, BH growth is 
self-regulated, resulting in a tight scaling relation between 
BH mass and spheroid mass. Following a merger, BH are 



allowed to grow until they reach a critical mass, at which 
the luminosity emitted by the AGN is sufficient to power a 
pressure-driven outflow that is able to slow down subsequent 
accretion, and eventually to unbind the gas in the galaxy. In 
this picture, BH accrete at the Eddington limit until they 
reach this critical mass, after which, in the "blowout" phase, 
the accretion rate declines as a power law function of time. 
The FID model reproduces the observed number densities 
of AGN at z ~ 2 over a wide range of bolometric luminosity, 
but in its original form it does not reproduce the observed 
downsizing trend: at low redshift, the FID model overpro- 
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duces the number of bright AGN and under-predicts the 
number of moderately luminous AGN, while at high red- 
shift, this trend is reversed. Therefore, we extended the FID 
model by considering different modifications to the physical 
recipes for black hole growth, and investigated their effect 
on the AGN evolution in different luminosity bins. We sum- 
marize our main findings: 

1. A sub-Eddington limit dependent on the cold 
gas fraction of the host galaxy (VE model): The 

FID model overproduces luminous AGN at low redshift 
(z < 1). We found that introducing a sub-Eddington cap 
on the BH accretion rate which was dependent on the gas 
fraction of the progenitor galaxies resulted in improved 
agreement with the observations in this regime. In our 
model, the gas content of massive galaxies is strongly de- 
pendent on the Radio Mode AGN feedback, which sup- 
presses cooling in massive halos at low redshift. Low cold 
gas fractions may retard the loss of angular momentum 
due to smaller viscosity and thus, the cold gas flow onto 
the central black hole may be suppressed. However, we 
found that we only obtained a good match to observa- 
tions when we implemented this accretion rate cap at 
z ^ 1. Extending the same gas- fraction dependent limit 
to all redshifts results in too few of the brightest AGN 
at redshifts between 1 < z < 4. One possible explana- 
tion might be that the evolution of the cold gas content 
in the FID SAM may not be correctly modeled. This 
could be due to incorre c t or simp l ified cooling recipes 
iHirschmann et ail (|20 1 if ) ; iLu et all (|201ll ) or inadequate 
star formation recipes (Caviglia & Somerville, in prep). 

2. An additional BH accretion triggering mecha- 
nism due to disk instabilities (DI model): The sec- 
ond main result of our study is that gas accretion onto 
black holes due to disk instabilities seems to be a non- 
negligible trigger mechanism for moderately luminous 
AGN with black hole masses 10 7 M© < M. < 1O 8 M 
at low redshift. This increases the number of AGN with 
bolometric luminosities 10 43 erg/s < Lboi < 10 45 erg/s, 
i.e. Seyfert galaxies, by about one order of magnitude. 
Our results therefore favor a 'hybrid' picture in which ma- 
jor merger events are the main driver of luminous AGN, 
especially at high redshift, while disk instabilities are the 
main mechanism powering moderately luminous Seyfert 
galaxies at lo w redshift (consistent with t he picture sug- 
gested by e.g. iHopk ms fc Hernauistll2009l ). Note that the 
studies based on the Munich and the Galform model 
have come to rather different conclusions (see discussion 
below). 

3. Heavy seed black holes with a halo mass thresh- 
old for seed formation (SH model): Our FID model 
does not produce enough very luminous QSOs at high 
redshift (z > 5). Our third, main conclusion is that a 
heavy seeding mechanism for black holes provides a phys- 
ically motivated way to increase the predicted number 
density of very luminous QSOs at high redshift, in line 
with observations. This is in ag reement with th e stud- 
ies of IVolonteri et al.1 (|2008T l and I Volonterll (|2010T ). which 
showed that either massive black hole seeds are required 
or less massive seeds have to accrete at super-Eddington 
rates. Unfortunately, current observational constraints 
are not sufficient to favor one of these possibilities. This 



may be possible with the next generation of planned X- 
ray missions (e.g. WFXT, IXO). As an implicit conse- 
quence of massive black hole seeds, we have also assumed 
that seeds are not able to form in dark matter halos be- 
low 2 x 10 Mq, as cold gas might not be able to col- 
lapse in these halos d ue to their shallow potential wells 
l|Volonteri et al.l 120081 ). This simultaneously reduces the 
number of faint AGN at high redshift, in better agree- 
ment with observations. 

We find that the FID model with a combination of the 
above modifications (=VEDISH model) can reproduce the 
downsizing trend fairly well and is e.g. able to m atch the 
observational compilation of iHopkins et alj (120071 ) . More- 
over, we have shown that the additional modifications do 
not change basic galaxy and black hole properties at z — 
significantly and thus, our model galaxies are still in agree- 
ment with the observed local stellar mass function, black 
hole mass function and the black hole-bulge mass relation. 
However, as in other semi-analytic models we find that at 
high redshift (z Si 1) the number of low-mass galaxies at 
high redshift is overestimated and the number of high-mass 
galaxies may be underestimated. These predictions are in- 
sensitive to our modifications to the black hole seeding and 
growth recipes, but of course our predictions for AGN num- 
ber densities are impacted by these discrepancies in correctly 
predicting the galaxy population. 

For the "best-fit" VEDISH model we find that the 
peaks of the distributions of Eddington ratios move to- 
wards smaller values with decreasing redshi ft in quali- 
tativ e agreement with obse r vational studies dVestergaardl 
20031: iKollmeier et al.1 120061: iNetzer fc Trakhtenbrotl 120071 ; 



Kelly et al l l2010l : ISchulze fc Wisotzkil l20ld ). Additionally, 
the results of the VEDISH model are in excellent agreement 
wi th the observed accretion limit in the Lboi — Mm -plane 
of ISteinhardt fc Elvis! (|20ld ). In our model, this behavior 
results from the gas-fraction dependent accretion rate cap 
that we applied. 

Despite the success of the final VEDISH model in re- 
producing the bolometric AGN luminosity function from the 
Hopkins compilation, the number of faint AGN still appears 
to be overestimated at high redshifts 4 > z > 2. One possible 
explanation for this discrepancy is a r edshift-dependent dus t 
obscuration, which was neglected in IHopkins et all l|2007h . 
Therefore, we have compared our model results directly to 
the observed AGN LF in the soft- and hard X-ray ban d by 
assuming the obscuration model from Ha singeJ (|2008h for 
the soft- and also for the hard X-ray luminosities, as even 
2 — 10 keV X-ray surveys will miss a significant fraction 
of mo derately obscured AGN and nearly all Com pton-thick 
AGN dTreister et al.ll2004l ; [Ballantvne et alii 20061 ). With our 
VEDISH model including obscuration effects we achieve ex- 
cellent agreement with the observed AGN LF in the soft and 
hard X-ray bands for the whole luminosity range at all red- 
shifts. Therefore, our results suggest that the obscured and 
Compton-thick AGN missed in deep X-ray surveys likely 
constitute a significant fraction of the total AGN popula- 
tion at all luminosities, in particular of moderately lumi- 
nous AGN at high redshift. However, we note that although 
obscuration effects contribute to the late peak in the num- 
ber densities of moderately luminous AGN, in our model 
they cannot account for the observed downsizing trend com- 
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pletely (see FID model and dust obscuration in Figs. [9] and 
E}. 

Comparing our conclusions with those of some previous 
studies (see Section [ 2l), we summarize the following main 
points: iMarulli et al.l ( 20081 ) found that a model based on the 
original Munich semi-analytic code and BH growth recipe 
l|Kauffmann fc Haehneftl l200d ; ICrotonl [ioOrj ), and with a 
light curve model similar to the one adopted here, under- 
produced luminous AGN at high redshift and slightly over- 
produced them at low redshift [z ~ 0), and underproduced 
low-luminosity AGN at low redshift (z ^ 1). As a result, 
they introduced an ad-hoc scaling of the black hole accre- 
tion efficiency with an explicit redshift dependence. With 
this "best fit" model, they were able to reproduce the bolo- 
metric AGN luminosity function fairly well except at very 
high redshift [z > 4). Thus, their results are qualitatively 
very consistent with ours — in effect, they found that the 
same qualitative modifications to the model were necessary, 
and achieved this via an explicit modification of the accre- 
tion efficiency, while we have tried to achieve this by adding 
more physically motivated effects. However, contrary to our 
study, they find that they do not need AGN activity trig- 
gered by disk instabilities, but can explain all AGN activity 
with a merger model. 

In contra st, in the study based on the Durham Gal- 
form model l|Fanidakis et al.ll2010T l. disk instabilities are 
the major triggering mechanism for black hole activity for 
all luminosities at all redshifts. With their model they can 
reproduce the AGN luminosity function at all redshifts: 
the low-luminosity end of the AGN luminosity function is 
mainly due to their ADAF model (cold gas accretion onto 
the black hole in a hot halo), whereas in our model it is 
due to accretion in the power-law decline regime as well 
as black hole accretion from disk instabilities. At high red- 
shift, they predict that super-Eddington accretion is respon- 
sible for producing luminous AGN. At this point, it is not 
clear whether a heavy seeding mechanism (as assumed in 
this study) or super-Eddington accretion describes the cor- 
rect physical processes. For that question, the improvement 
of high-redshift AGN observations will be of crucial impor- 
tance, both by enlarging the current high-z AGN samples 
and by reducing the current uncertainty originating from 
incompleteness problems. They mainly attribute the down- 
sizing trend to dust obscuration effects. 

Thus the main triggering mechanism of most of the 
black hole growth in the universe remains a major open 
question. It is clear that with the allowed freedoms from 
various modeling and observational uncertainties, it is pos- 
sible to reproduce the number densities of AGN (luminosity 
functions) as a function of redshift within a range of scenar- 
ios, from a pure merger scenario, to a hybrid merger+disk 
instability scenario such as the one we have suggested, to 
a pure disk-instability driven scenario. Hopefully, this issue 
will be clarified by further studies of the morphology of AGN 
hosts, now possible out to the peak of luminous AGN activ- 
ity z ~ 2 with Wide Field Camera 3 (WFC3) o n Hubble (e.g. 
Schawinski et al. I l201ll : iKocevski etafl 12012] ; H osario et al.l 
201 If ), and the enhancement of AGN activity in close pairs 
over a range of redshifts and luminosities. We intend to make 
a more quantitative comparison of our predicted host prop- 
erties with observations in a future work. 



Overall, we conclude that the models presented here 
provide a plausible attempt to understand the complex sce- 
nario of black hole and galaxy co-evolution, and to predict 
the downsizing trend within the framework of hierarchical 
clustering. However, there still remain many uncertainties 
in modelling the formation and evolution of black holes, 
which hopefully ongoing and future observational facilities 
will help to constrain. 
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